INTEGRAL APOLLONIAN CIRCLE 
PACKINGS AND PRIME CURVATURES 

J. BOURGAIN^ 



Summary: It is shown that any primitive integral ApoUonian circle packing captures 
a fraction of the prime numbers. Basically the method consists in applying the cir- 
cle method, considering the curvatures produced by a well-chosen family of binary 
quadratic forms. 



Introduction 

In this paper, we pursue a line of research initiated in [GLMWY] and [S] on the 
arithmetical properties of integral Apollonian circle packings (ACP for short) in the 
plane. The reader is also referred to [B-Fl] for certain background material. 

Throughout the paper, we consider bounded ACP's which are primitive, meaning 
that all curvatures of the circles in the packing do not share a factor greater than 
one. Let us recall that the set of curvatures in a given packing P is obtained by 
action of the Apollonian group A on the root quadruple (a, 6, c, d) of co-prime integers 
a<0<b<c<d, a + b + c>d. The group A is a subgroup of the orthogonal group 
associated to the Descartes quadratic form 

Q{xi,X2, X3, X4) = 2{x\ + X2 + X3 X4) - {xi -I- a; 2 + a; 3 -I- x^)'^ 

whose vanishing is tantamount with xi, 0:2, 0:3, 0:4 being curvatures of mutually tangent 
circles. The group A is generated by the matrices 
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The basic problem is to describe the set of curvatures appearing in a given packing 
P; the ultimate hope is to establish a local to global principle cf. [F-S]. More modestly, 
[GLMWY] put forward the 'positive density' conjecture, according to which the set 
of curvatures in an ACP form a subset of Z of positive density. Following up on 
a technique proposed by P. Sarnak, this problem was solved afl&rmatively in [B-Fl] 
(a slightly stronger result is obtained in this paper; see Theorem 1 and the Remark 
following its proof). Using a result due to Iwaniec on representing shifted primes by 
binary quadratic forms, Sarnak also pointed out that any primitive ACP produces at 
least CjY^^Y)372 distinct prime curvatures at most X, for X ^ cxd. Based on new results 
on the representation of integers by binary quadratic forms of large discriminant, 
previous lower bound for the number of prime curvatures is improved further in 
[B-F2] to at least t^i^I2—- 

The main result in this paper gives the correct order of magnitude. 



Theorem 2. Given an integral primitive ACP, there is a positive c, such that for X 
large the number of prime numbers less than X which are curvatures of circles in the 
ACP, is at least Cj^j^ (with c> an absolute constant). 

Compared with the arguments due to Sarnak and refined in [B-F2], that are based 
on Iwaniec' theorem and representations by individual quadratic forms, the strategy 
used here is different. Our approach consists in introducing a generating function by 
considering the collected contribution of suitable families of binary quadratic forms 
(constructed in §0, §1 of the paper). These generating functions can then be analyzed 
using the circle method (in a rather standard way), to the extent of providing a main 
(arithmetical) contribution with an error term. In particular, we are able to establish 
Theorem 2 (relying also on the so-called 'majorant property' for the set of the prime 
numbers, established in [B], [G]). The technique applied here may be organized better 
as to allow a treatment of the major- arcs contribution by spectral methods (using the 
spectral analysis for the full Apollonian group), in the spirit of [B-K]. This leads to 
better error terms and statements that come close to a local to global principle. That 
program is pursued in the forthcoming paper [B-K2]. Let us also mention the paper 
[F-S] that gives evidence for the only congruence obstructions to appear (mod 24). 

In this discussion, we should cite the paper [K-0] , where counting results for the cur- 
vatures, with multiplicity, are obtained based on spectral techniques (see also [BGS]). 
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In particular, it is shown in [K-0] that in any ACP P the number of curvatures at 
most X is of the order 

for X ->cx) 

with 5 = 1,30068.. is independent of the packing. This amounts also to the number 
of quadruples bounded by X in the orbit of the root quadruple under the Apollonian 
group A. 

Let us briefly recall how binary quadratic forms enter the analysis (see [S]). While 
^ is a 'thin' (non-arithmetic) group, its subgroup Ai — {82,83,34) (= stabilizer of 
xi) and similarly ^2)^3)^4 are arithmetic. More precisely, considering the map 

y = {y2, 2/3, 2/4) = {x2, xz, xa) + (a, a, a) 

the afiine action of Ai on (x2,X3,X4) is conjugated to the action of a finite index 
subgroup r of OgCZ), g denoting the quadratic form 

giv) = y2 + yl + yl- '^y^ys - 22/22/4 - 22/32/4- 

By a further coordinate change 

A = y2,B= ^{yi - ^3 + 2/4), = ^4 

transforming g{y) in the quadratic form A{A, B,C) — — AC, V is conjugated to 
the subgroup of Oa(Z) generated by the reflections 
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Consider the spin double cover of 5'Oa(Zi) realized as image of GL2(Z) under the 
homomorphism 



P ■ 



a 13 
7 5 



1 



I-)- 



a5 — 



a(3 ad + ■yd 



2PS 


s\ 




"1 -2" 
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and hence the 



with kernel ±7. Then p-^(50a(Z) n F) contains 

principal congruence subgroup A (2) of 5'L2(Z). 

It turns out that the set of values of 2/2 = A, ys = A + C — 2B,y4 = C contains 
at least those of A with {A, B, C) ranging in an orbit p{SL2{'L)){Aq,Bq,Cq)\ i.e. the 
integers represented primitively by the binary quadratic form 



Aoa^ + 2^00:7 + Co7^ with (a, 7) = 1. 
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The preceding provides an explicit recipe to produce curvatures in a given packing P. 
Assume (oq, cq, do) e 5 = S{P) = A{a^ b, c, d)* and set 

Ao = ao + bo, 2Bq = oq + 6o - cq + rfo, Co = ao + cio- 

Then aU integers represented by the quadratic form 

Aqx^ + 2Boxy + Coy^ - ao with x,y eZ, (x, y) = 1 

appears as curvatures of circles in the packing P. 

This observation made in [S] plays a key role in [B-Fl] and also in the construction 
of an appropriate family of binary quadratic forms described in §0, §1 of this paper. 

(0). Preliminary construction of a set of curvature quadruples 

Let -Ri be a large integer and denote Sr^ the set of quadruples (a, b,c,d) E S — S{P) 
of the Apollonian packing P satisfying 



Thus 



Let 



max(|a|, \b\, |c|, \d\) ~ Ri. 

13 



\Sr^ \ > Rl with 5 > 



10 



(0.1) 
(0.2) 



i?2 = i?f. (0.3) 

Given (a, b, c, d) e Sr^ , let A = a + b, C = a + d, 2B = a + b — c + d and consider the 
set of integers 

Sa,b,c,d = {Ax'^ + 2Bxy + Cy^ — a; x,y E Z,0 < x,y < R2 and {x, y) = 1}. (0.4) 
Recall that (A, S, C) = 1 and = AC — B^ (by Descartes' equation). 

As explained in [B-Fl] the set Sa,b,c,d is contained in the set of curvatures produced 
in the orbit of (a, b, c, d) under group elements of Ai — {S2, S^, S4) of norm bounded 

/"\ 

by Denote S{a, b, c, d) C Sr_^r2 n{S2, S3, S4) 



c 

\dJ 



a set of quadruples {a,b' , c' , d') 



in one-to-one correspondence with Sa,b,c,d by projection on the ^'-coordinate. 

Thus for each < x,y < R2, {x, y) = 1 there is some gx,y € ('S'2, Ss, S4) such that 



b 



c 

\dj 



S(a,b,c,d) C ^gxy ^ ■,0<x,y<R2 and (a;,y) = l|. 



(0.5) 



Obviously, if we fix a;, y, all quadruples g^y 



b 
c 



are distinct and hence 



l5(a,6,c,d) < R 

(a,b,c,d)&SR^ 



(0.6) 



Given b', it follows from Descartes' equation that 

7r^\b')nSj,^j,.<^Rl+'Rl 

and (0.6), (0.7) imply that 

{a,b,c,d)£Sii.^ 



(0.7) 



(0.8) 



To each (a, b, c, d) e Sr^ , associate the distribution Xa,b,c,d on Z obtained as image 
measure of 

[0 < x,y < R2:,{x,y) = l,fa{x,y) ^ RiRl and {fa{x,y), p) = 1] 



P<Rp 



under the map 



(a;, y) ^ /„(a;, y) = Ax^ + 2Bxy + C?/^ - a. 



Hence, suppAa,6,c,d e Sa,b,c,d, ||Aa,6,c,d||oo < -Rf and elementary sieving implies that 
certainly 



^ ||Aa,6,c,di|l ^ -^2- 



(logi?2)2 

Also, from sieving, we obtain that for all q e Z+ 



(0.9) 



z=u(mod q) 



E >'a,b,c,di^) ^ (g + ll'^a-''.c,d||l 



(0.10) 



(in the argument, we distinguish the cases q > (logi?2)^^^ and q < (log i?2)^°'^; we 
only need a crude estimate for our purpose) . 



Define 



A = E K,b,c,d 

{a,b,c,d)€SR^ 
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(0.11) 



which is a distribution on [h' e b' RiRj]- From (0.8), (0.9) 

'^'(logi?2)^ 



Sr,\ ,^ < ||A||i < \SrA-RI and ||A||oo « Rl+'Ri (0.12) 



By construction, z G suppA has no prime factors less than 7^2° and from (0.10) 
obviously 

E >^i^) ^ (- + (0.13) 

z=u(moa q) 

Let r]a,b,c,d be a distribution on 5(a, 6, c, d) which image measure under projection on 
the ^'-coordinate equals \a,b,c,d and set 

V= Yl ^«Ac,d. (0.14) 

{a,b,c,d)eSRj^ 

Hence 

\\v\\i = \\M\i- (0-15) 

Since clearly ||r7a,6,c,d||oo < ||Aa,6,c,d||oo < Ri, it follows from (0.6) that 

||77||oo < (0.16) 

(0.13) may be rephrased as 

E Via',b',c',d') < + "ll^lli for aU q e Z. (0.17) 

6'=tt(mod q) 

Next, we replace the distribution r] hj a subset C C iSjj^^2 which we construct 
probabilistically by selecting {a' , b' , c' , d') e C with probability 

5r]{a', b', c', d') < 1 where 5 = R^^ (cf. (0.16)). 

By (0.12), (0.15), we obtain 

\C\ ^ 5.\\r]\\i > R^^ISr^I (0.18) 

and also 

E ~ ^ll^lli + ^'ll^lli 11^11°° < (0.19) 

6' 
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with Cb' denoting the fibers of C. 

From (0.17) and standard large deviation inequahties, we deduce that for g e Z_|_, 

\{{a', b', c', d!) eC;b' = w(mod q)}\ < 

6 Yl v{a',b',c',d') + cV^^i{SMi)-^ < 

6'=u(mod q) 

^ + ^T\C\. (0.20) 



q i?2 

Relabehng the quadruples, we obtain a subset C C Sii,R= R1R2 with the following 
properties 

(0.21) \C\ > R^-To > i?6/5 

(0.22) o ~ and |c|, \d\ < R for (a, 6, c,d)eC 
(0.23) EalCaP<i?-^/'|C|2 

(0.24) For (a, b, c, d) e C, a has no prime factors less than R'^ 

(0.25) \{{a,b,c,d) eC;a = u{mod q)}\ < {q-^ + R-^)3^\C\ for all q e Z+ 

(1). Introducing a family of quadratic forms 

Let C C Sr be the set constructed in §0. Let A — TTa{C) C 
To each (a, 6, c, d) e C, we associate again the binary form 

f{x,y) ^ Ax^ + 2Bxy + Cy^ (1.1) 

with 

A = a + b,C = a + d,2B = a + b-c+d, disc/ = -4a^ = - S^. (1.2) 

Thus iA,B,C) = 1. Since |C| < R and \A.C\ > i?^, it follows that 

|A|,|C| ^R. 

Denote J-" the family of quadratic forms (1.1) obtained from C and by J^a C J-' those 
obtained from Ca- We show that if we fix the discriminant, the number of equivalent 
forms in is 0{1). 
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Thus if Ax'^ + 2Bxy + Cy'^ and A^x'^ + 2Bixy + Ci y^ are equivalent, then 



Ai = a^A + 2a-fB + 7^(7 
Si = a^A + (a5 + p-i)B + 7^(7 



for some e 5L2(Z). 



Thus 



Hence 



a 

and since \ Ai\ ~ i?, 



7^ < ^ < 0(1) 



« + |7 



<0(1),|«|<0(1) + 0(1)^<0(1). 



Similarly, we see that, < 0(1). This shows that at most 0(1) quadratic 

forms obtained from Ca are equivalent. This proves our claim. 

As a consequence, we obtain that for fixed a and M 

#{{f,h,x,y,x^,y,) e:FaXTaX [l,M]^;f{x,y) = h{xuyi)} « {RMyM^\J^a\. 

(1.3) 

Indeed, f E J^a and x, y. The integer z — f{x, y) is at most i?M^ and is represented 
by at most 2'^^^^ classes with discriminant —Aa?. From the preceding, there are at most 
0(1) 2'^^^) possibilities for /i e IFa and for each fi, fi{xi,yi) = z holds for at most 
[RMy values of {xi,yi). This establishes (1.3). 

2. Application of the circle method 

Let R and be as in §1. Let P be large and assume R < (logP)*-^. 
For / e J", let a;/ be the image measure on Z of [1, P]^ under the map 

(a;, y) H> fa{x, y) with f^ = f-a, -4a^ = disc / 

under restriction {x,y) = 1. 

Actually it is technically more convenient to consider the image measure of 
^{x y)=i'l{^) ®7(p); where < 7 < 1 is a smooth bumpfunction supported on 
[0,1]. 
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Note that suppw/ C [— i?P^, RP^] (by construction of JF) and we may assume 
suppw/ C [0,i?P^] (by assuming ^ > 0). Define 



= where = Yt\^ ' 



The co-primahty condition {x,y) = 1 leads to technical complications (with no 
effect on the basic scheme of the argument). 

Fix some integer B = (logP)^° and replace the restriction l(x,y)=i by 

d\{x,y),d<B 

This expression equals 1 if (x^y) = 1, vanishes if 1 < (x^y) < B and is bounded by 
T((a;, y)) . Hence this replacement introduces an error at most T(^{x.y))l(^x,y)> b <^ 

^(logP)^ for the counting functions Uf and u) (in £-^(Z)-norm), which is harmless. 
With the above modification, we obtain 

u;i9)=J2f^id){T E 7(|)7(|)e(/.(a.,y)^)} 

T d\{x,y) 



d<B 



and 



\0J 



E-^'^'lE E 7(^)7(|)e(/.(x,,)9)f. (2.1) 



d<B JF d\(x,y) 

Fixing d square free, we have to analyze the expressions 

SUO) = Y. E ^{^)^{^y{iAx' + 2Bxy + Cy'-a)e). 

T d\{x,y) 

For q < RP, denote for (g, b) = 1 



(2.2) 



M{q,b) 



e- 



< 



qRP 



C T 



Our main concern is to obtain suitable bounds on 



E E 



\s^{e)\^de. 
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(2.3) 



Lete=^ + ifeM{q,b),\ip\<-^. Then 



faix,y)9 = faix,y)^ + faix,y)(p 



and 

E 7(1)7(1) e(/a(a;,y)^) = 



d\{x,y) 



E iuk,i)^)[ E </«(-^i/H7(|)7(|) 



o<fc,^<9 ^ X = k,y = £(mod 

where do = {d^q),d = d^di. 
We distinguish 2 cases. 

Case I. q< P. 

Rewrite the second factor in (2.4) as 

'rq + k\ fsq + 



r s<^ 
rg+fc=sq+^=0(mod d\) 



(dropping a multiphcative factor) 

■r'diq + k'\ /s'diq + £'' 



E 7(^^^)7(^^)e(/(r'd., + .',.'d., + 0^) 



with < k', £' < diq; k' = k{mod q),e' = £(mod q), k' = £' = O(mod di). 
Prom the Poisson summation formula, we obtain 



— J,2 E Jf{Qdi,m,n,(fi)eq{—mdik — ndi£) 
5 "1 ^ ^c'y 



m,n€Z 

where di = 1 (mod q) and 



Jf{q,m,n,(fi) = JJ l{^^l(^^^e{f{y, z)(fi)e(J^y + ^z^dydz. 
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Note that, by stationary phase 

1"^^^^'™'"'^)" - 1^1 |disc!(/)|V2 ^ -i^y (2.7) 

Also, in (2.5) the significant contributions come from values m, n satisfying 

MM<|V/|.M<B.P.^ (2.8) 

hence |m|, \n\ = 0{di). 

From (2.7), there is an obvious bound on (2.5) and on (2.4), 

(2.2). Substituting in (2.3), we see that the contribution of \(p\ > (p^ is at most 

< —{logP)-^ provided >(logP)^^. 

1 5 . 

If Q < Pto , we may therefore restrict \(p\ < P~'s and obtain m = n = in (2.8). Thus 
we distinguish the ranges 

Q < Pto . Then (2.4) is replaced by 

Y: e[Uk,£)^)]jf{^) (2.9) 



0<fc/<g 



contributing in (2.3) for 



1 \\2 



do\{k,£) 



(where |c/| < 1). 

Summing over d, we obtain 

p2 



(2.10) 



^E^o+lE E |E^/^/.'^o(?,^o,o)|'} (2.11) 

do q~Q (b,g) = l .F 



do\q 



denoting 



Sf{q,b,m,n) = \ J] eq{hfa{k,t)-mk-n£) (2.12) 
5/,d„(g,6,m,n) = 4 E eq{bfa{kj)-mk-n£) (2.12') 



9^ 

0<fc,£<g 
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and where do < B is a, square-free divisor of q. 

Pto <Q<p 



By (2.5), (2.4) becomes 

i XI Sf^do{Q,b,mdi,ndi) Jf{qdi,m,n,(p) (2.13) 



^ |m|,|n|<0(di) 

and the contribution in (2.3) is at most 

p2 



^'^V.^^^ E{ E \j2'^f^i^f^do{Q,b,md^,nd^)\] (2.14) 

H do,di,m,,n ^ — ' L ^ — ' I ^ — ' I J 

g~Q (6,g)=l ^ 
do\q 

with |c/,q| < l,d,di square free, {do,di) = l,dodi < B and m,n = 0{di). 

Case II. P <Q <PR 

Since for ^ = ^ + G A^(g,^) 

/„(a;,2/)^ = /a(a;,y)^ + /„(x,2/)^ with |/„(a;, y) | |(^| < i^.P^.-^ < 0(1) 
we may replace ^ by | (dropping (p) and Su>id) becomes 



d\{x,e) 



Thus the inner sum in (2.15) equals 

E ^,(f/«(<ii*,<iiQ)7(^)7(^) 

do\{k,£) 

and completing the sum, we obtain 

eq{-ab). ^ eq{djbf{x,y) + ux + vy) 

0<x,y<q 
do\ix,y) 



\u\,\v\<^di 0<x,y<q 



= q^ ^ Sf^doiQ^b, diu, div). (2.16) 
|u|>|<Sdi 
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The contribution to (2.3) may be bounded by 

E E I E J2^f^doiQ,b,d^u,d,v) (2.17) 

qr^Q {b,q)^l u,v<%di ^ 

with d^^dx square free, ((io,c^i) = 1 and d^dx < B. 

(3). Evaluation of the Gauss sum 

Analyzing further (2.11), (2.14), (2.17), we evaluate expressions of the form 

E Sf{q,b,u,v)Sf^{q,b,Ui,Vi) (3.1) 

and also 

E ^f,doi(lib,u,v)Sf,,do{Q,b,ui,vi). (3.1') 

ib,q) = l 

Consider first (3.1) 

Factoring q as a product of prime power p'^, (3.1) factors correspondingly. Recall 
that ^ 

Sf{q,b,u,v) = ^ E eq{bf{x,y) + ux + vy-ba) (3.2) 

^ 0<x,y<q 

with /(x, y) = Ax'^ + 2Bxy + Cy^, {A, B, C) = 1. 

Let g = p'' be a prime power. We may assume {A,p) = 1. 

Write (assuming p 2; for p — 2 there are some extra technicalities that we omit 
here and to which we will return in §7) and using the notation ~ for the multiplicative 
inverse (mod q) 

bf{x, y) + ux + vy — ab = b{Ax^ + 2Bxy + Cy^) + ux + vy — ab 
= bA{x + BAyf + ba'^Ay'^ + ux + vy - ab 

= bA{x + BAy + 2bAuf + ba'^Ay'^ + (v - BAu)y -WAu^ - ab (mod q). 
Hence, by Gauss sum evaluation, we obtain (cf. [BEW]) 



^/(g,6,w,^;) ~ )e,(-4A6w2 -a6)[ ^ e^iba" Ay"^ + {y - BAu)y)^. (3.3) 
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Let ^ = I with (a, q) = 1. Thus writing y = y + zq,0 < y < q, < z < i = {q,a^) 



0<y<q 0<y<q 0<z<{a'^,q) ^ ^ V ) 

and (3.4) = unless (a^, q)\v — BAu, in which case we set 

V - BAu = {v- BAuy.ia^, q) 

Thus (3.4) becomes 

(a^g). eg{bAdy^ + {v-BAuyy) 

0<y<q 

= (a^?) Yl eq{bAa(y+{Av- Bu)'^2Mf -mA{{Av- Bu)^f) 



(3.4) 



0<y<q 



{a\ q){qy/' " ^ ((^^ " BuD') ■ 



(3.5) 



Hence, from (3.3), (3.5) 

(a2,?)V2 ^bA./bAh 



Sf{q,b,u,v) 



(^^) (^£|^^e,(-OT«2 _ a6)e,~(-46SA ((At; - Sw)~)') 
eg(-4]45M2 - a6)eq(-45a]4 ((Av - 5w)~)^) if ^ 1 



(a2,g) 



= < 



(a2,g)^ 2 (3.6) 

-^^ — — jeq{- AAbu^ - ab) iiq\a^r ' 



It follows that 

^3_;L) ^ (a^g)^(al,g)^ 



Y eg((4Aiu^ -4Au2)6 + (ai -a)6^E(6)Ei(6) (3.7) 



(b,Q) = l 



where £■(&) is either 



(t) 



or Cgf — 4yla ((Av — S'u)^) 6) with g = , „, 
g / ^ ^ ^ ^ (g, a^) 
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(3.8) 



where we assume {q,a'^)\Av — Bu and similarly for Ei{b). 
Thus the sum in (3.7) is of Kloosterman or Salie-type. 
There is the following elementary estimate (which suffices for our needs) 

I Yl epr{cb + db)\<ip'-f/\p',c,d)^/^ (3.9) 

{b,p)=i 

0<6<p" 



and also 



J2 {^)eAcb + db)\<{p^f/\p^,c,d)'/\ (3.10) 



{b,p)=i 

0<f.<p" 

Hence we may state the following bound for q — p^ 

(3.7) < {a^qy^alq)Hq,a-a^f/\-'>/^ (3.11) 

and also, for a = ai 
(3.7) < 

q-^'^{a^,q){q,lAiul-^u^ + {q, a^)M^{{AiVi - Bimyf - {q,a^)UE{{Av - BuDY^^ 

(3.12) 

where {a^ , q)\Av - Bu, {a^,q)\AiVi - Bmi and a = ^ , {Av - Bu)-" = . 
Consider next (3.1') for which there is again factorization. 
If q = p^ and p\do, then 

Sf,p{p'',b,u,v) ^ p~'^''epr-{-ba) ^ epr-i{bf{x,y)p + ux + vy). (3.13) 

Note that by assumption (0.24) and since p < (logP)^°, {a,p) = 1. 

We distinguish several cases 
(3.14) r = 1 

Then |(3.13)| = p~^ and |(3.1')| — {p — l)p~'^, which is certainly bounded by 

p-i (3.11) (3.15) 
15 



and 

p-i (3.12). (3.16) 

(3.17) r > 2 

Clearly (3.13) vanishes unless p\{u, v). Writing u = pui,v = pvi, 

(3.13) = epr{-ba)p-^p-^^^-^^ [ ^ epr-.{hf{x, y) + mx + v^y)\ . (3.18) 

0<a;,j/<p'"-2 

If r > 2, repeating the analysis of the exponential sum with q replaced by qi = p~'^q, 
we obtain instead of (3.6) 

1 



-eg, ( - 4Abu'^ - 4ba'^A{{Av - Buy) )eq{-ah). (3.19) 
This gives for (3.1) the bound 

(?,a-ai)VV'/^ (3.20) 

Factor 

^= n WpWp"^ 

(p,do) = l p\do p\cLo 
r>l 

and consider the corresponding factorization of (3.1'). Apply (3.11), (3.12) if (d^p) = 1, 
(3.15), (3.16) if p\do, r = 1 and (3.20) if p\dQ, r > 1. 

(4). Estimation of (2.11), (2.14) 

Expressing the square of the inner sum and carrying out the summation over 6, we 
evaluate (3.1), (3.1') by (3.11), (3.15), (3.20). 

Hence ^or Yli{h q)=i\" '\'^ '^^ obtain the bound 

Q-'" {a\qf'\alqf'\q,a-arf'\\{p-y'). (4.1) 

p\do 
r=l 

Substitution in (2.11) gives 

sji S EE('''.9)*«9)*('>-''-9)i(np-'^')(nf'^') 

fJl&J^ do qr^Q p\do p\do 

^^r^ E Y.(^',Qi)kal,qi)ka-a^,q,)Ul;' (4.2) 
do\q 
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where q = qido2 and cio 2 = H P- Note that do\qi and hence for given qi, there are 

' ' p\do 

r>l 

at most min(gf , i?^"*"^) possibihties for (io,i,<^o,2 and q. 

Specify qi ~ Qi where Qi < Q < m.m{BQi,Ql) and replace (4.2) by 

p2 



miniB^Ql) ^ ^ (a^, gi)^ (a?, (a - ai, (4.3) 



Next, we evaluate 

p2 



The contribution for a = ai may be bounded by 

p2 1 



^ ' ' oG^ ' ' a (4.5) 

For a 7^ a 1, we obtain 

p2 



L E l-^a||-F„J5^(a^g)V2(a?,g)V2(a-al,g)V^ (4.6) 



Clearly, for fixed a 7^ ai 



^(a^,)V2(a?,,)V^(a-ax,g)V4< E ^'^'^^^'1^ ^ (^.7) 

[d,di]<Q 

since a, ai are pseudo- prime by (0.24). 
Also 

J2 (a^ qy/\al qf'\a - a^q)'/' < E («' 9)(a - ai, g)^/^ 

q~Q q^Q 

d|o,di|oi,d |o— ai V / 

[d,di,d']<Q 

17 



Denote / = (d^di). Thus f\d' and since / ^d,f ^di,f ^d' are pairwise coprime, 
[d,di,d'] > f~^ddid' and therefore 

(4.8) <Q.(a,ai)3 J] {d')-'/\ (4.9) 

d'\a—ai 
d'<Q 

From (4.7), (4.9) 

^"^■^^^TJF^ ^ l^a|.|-Fa,|min(l,g-^(a,ai)3 (i)^). (4.10) 

a,aieA d'\a-ai 

We distinguish two cases. 

Assume (o, ai) — A > Q^^ Estimate 

(4-io)<^^Ei-^«i- E E i-^-i- (4-11) 



aS-^ A|a,A>QiO"'' ai=a(mod A) 



Again, since a is pseudo-prime, A is restricted to 0(1) values, once a fixed. From 
(0.25), 

|{/ eT;a = u{mod A)}\ < A"^ \J^\ 



and hence 



(4.11) < (4.12) 



Assume (a, ai) < Q^^ . Then 

a,aiGA d'\a-ai 
d'<Q 



d'<Q 
p2 

< -B-<5"^- (4-13) 
-ft 



From (4.5), (4.12), (4.13), we obtain 



(4.4)<^(i?-H+Q-'°"\ (4.14) 
18 



We assume R ~ (logP)*-^ with C a sufficiently large constant. 
Since B ~ (logP)^", it follows from (4.14) that 



(4.3)<^(i?-Ug-^°"). (4.15) 
The same bound also holds for (2.14). 

(5). Estimation of (2.17) 

Expressing the square of the inner sums in (2.17) gives 



^1. 



RP 

qr^Q u,v,ui,vi<^di = ! 

do\q 

(5.1) 

where the inner sum is of type (3.1'). 
Note that here do, di are fixed. 

To estimate the contribution for a ^ ai, use again the bound (3.11) on (3.1) (which 
from previous discussion is always valid). We obtain 

f,heT q~Q 

< -Q'^RiB^+' < py^R^ < —R-\ (5.2) 
P R 

Next the a = ai contribution 

Let q = qiq2 with ^'2 = 11 factor (3.1') according to q = Yll,- If p\qi, the 

p\do ^ 
r>l 

bound (3.12) applies to (3.1') with q—p^. 

For p\q2, apply (3.11) which gives the bound g^^ on the g2-factor. Let gi ~ Qi, ^2 ~ 

Q2 (noting that, by construction, the number of g2-values is at most Ql)- Thus the 
contribution to (5.1) may be bounded by 

^^^^ E ^ E E E (8-3) (8.4) 

do,i\qi 
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with 

(5.3) = JJ Sf^doAP'^ib, diu, div)Sf,,do,i{p''ib, dim, divi) 

pWi {b,p) = l 

andc^o,l = ^p|do,r=l^'• 

Because {do,di) = 1, the factor di in (5.3) turns out to be irrelevant and we drop 
it for simphcity. 

Applying (3.12) for each prime p\q, we obtain 

p\qi 

{f jMul-'Um^ + (p^a^)4:47a((yll^;l - Si«i)~)^ - {f ,a^)^M{{Av - BuYf)"'^. 

(5.5) 

We distinguish several cases for the factors in (5.5). 
Assume 

(5.6) 

Multiplying with gives then 

lAu^a'^ + alAa(Av - Buf = JA^uja'^ + alA^{AiVi - Bmif (mod p[51+^) 
and since aa = l(mod j^r^) and Av — Bu = A\vi — Biui = (mod (a^,p^)) 

Au^a^ + A{Av - Buf = Axu\a^ + A^^A^vx - Bmif (mod ^[^1+^). (5.7) 
Hence, since = AC - B"^ ^ AiCi - Bf, 

f{v, -u) = /i(vi, -m) (mod plil+i). (5.8) 
Since 1^41, \B\, \C\, \u\, \v\ < R, it follows that either 

f{v,-u) = Mvx,-u) (5.9) 
or (5.8) can only hold for a set a of primes p\q such that 

l[p^i^+^<R\ (5.10) 

peer 
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Hence, if (5.9) fails, 



(5.5)<(a^9l)g7'/'Q;/^i^^ 

The contribution to (5.4) is at most 



(5.11) 



(5.12) 



Finally, the contribution of a = ai and {u,v,ui,vi) satisfying (5.9). We obtain 
B^+'QIqI^^ 1 ( P \^ 



R.P 



{qS[) E (5.14) 



with 



q^Qi aeA 



(5.13)= |(/,/i,'U,t;, wi,wi) e J"^ X -Jiu.v) ^ f{ui,vi)^ 

Recalling (1.3), it follows that 



and 



(6.13)<B'|Jb"|J-„| 



(5.14) < R'B^9!gl j^^j: 5:(a^,)|^.| 



QP (0.21) p 



(under proper assumption on R). 
Hence from (5.2) and (5.15) 



pii/io 



(5.15) 



P2 

{2.17)<—R--o. 



(5.16) 



(6). Minor arcs estimate 

From (4.15), (5.16), we obtain the following bound on (2.3) 

(2.3)<^(i?-ife+Q-^o"'). (6.1) 
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Take R = (logP)^°'° and 1< Qo < i?. 
Summing (6.1) over Q > Qo gives 



/ 



U U M(q,b) 
q>Qo (.b,q)=l 



Returning to the definition of Ai{q, b) in §2, we may further reduce the 
defining for q < QQ,{b,q) = 1 



M„{q,b) = 



q 



2 -, 



< 



RP'^ 



It follows from (2.4) and (2.7) that ioi 9 e M{q,b),9 = ^ + cp 



\SU0)\ 



< 



^R\(fi\ 



and therefore 



M{q,b)\Mo{q,b) 



^ ^ RQlq 



p 

which collected contribution is at most ^t^- . Therefore 



I 



\SM\'d9<^Q^''- 



n u u Mo{q,b) 

q<Q0 (b,q) = l 



where 1 <^ Qo < (log-P) is a parameter. 



(7). Contribution of the major arcs 

Let 

K Qoo < (logP)^ and Qi= Yl q. 

q\Qoo 

Let 

p2 
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and A a smooth even function on R, supp A C [—1, 1], < A < s.t. 

m 

Let V be the distribution on Z defined by 
Hence 



Estimate 



Q<Qoo 

where we apphed (6.7) with Qq replaced by Qqq and denote 



e- 



Qqo 
RP^ 



Moo{q, b) -- 
Since q\Qi and A is even 



'de 



(6.5) 
< 



max |1 — cos(mQiy?)|^ ^ d^p 



\m\<M 



< 



< 



6 

00 



qR^ 



and we obtain from the choice of R, Q 



00 



p2 

1 ^ ^ ,n-lO" 



1^ ~ 'S'(^,||2 < ^^QoO 
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Hence, by Parseval 



p2 

2 ^ /n-10" 



(7.7) 



Next, from definition of v and a;, clearly 



1 



i^(^) ^ \{{^^y) e [1,^^]^; = l,/a(2;,^) = n(mod Qi) and |/„(a;, < 

(7.8) 

Replacing the condition (x, y) = 1 by the weaker condition {x,y,p) — 1 for p < B 
introduces an error at most ^ with respect to the £^(Z)-norm. 

With this replacement, we obtain the following lower bound on (7.8) (since Qi,B <^ P^). 
XI |{^ - ^'^ ^ (^'2/'^oo) = 1 and fa{x,y) = n(mod (5i)| . (7.9) 



where (7.9) equals 



mes Isl, Itl < 



P 



mes 

Qoo<P<B 



f{s,t)- 



Isl, \t\ < 



n 



< 



P' 



pQi 



Qi QV 



f{s,t)- 



n 



< 



P' 



p'^Qi p'Qii 



(7.10) 
(7.10') 



From the assumption on f E J-', it follows that taking n ~ RP^ 



mes 



Isl, \t\ < 



f{s,t)- 



n 



< 



and similarly for the terms in (7.10'). Hence fiji) may be substituted by 
1 



(7.11) 



QiR 



^\{0 < x,y < Qi;f{x,y) = n + a(mod Qi) and {x,y,Qoo) = 1}\ 



R 



XI n [P '^\{^ ^^^y <P''A^^y^P) f{x,y) = n + a{m.oAp'')}\] 

(7.12) 



feT p<Qoo 

where Qi = Up<QooP' ■ 

Recalling (0.24), a E A has no prime factors less than i^sfe and hence (a, Qi) = 1 



Let 



f{x, y) = Ax^ + IBxy ^Cy^ eT where o? = AC - B"^ . 
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Fixing p < Qooj we may, since {A,B,C) — 1, assume {A,p) = 1 (the other cases are 
similar, except for p= 2, A,C even and B odd; we leave the adjustment to the reader). 

The equation 

fa{x,y) = n + a {mod p'') 

becomes 

{Ax + ByY + a^y^ = A{n + a) (mod p'') 

or equivalently 

+ = A{n + a) (mod p''). (7.13) 

We seek for a lower bound on the number of solutions of (7.13) with {x,y,p) = 1. 

Assume p > 2 (the local factor at p = 2 requires an additional congruence assump- 
tion on n at the place p = 2). 

Clearly, the number of solutions of (7.13) is at least 

p''~^\{{x,y) e {Z/pZy X Z/pZ;a;2 +2/2 = A{n + a){mod p)}\. (7.14) 
Using Gauss sums, we obtain 

|{(a;, y) e (Z/pZ)* x Z/pZ; x'^ + y'^ = A{n + a)(mod p)}\ 

= \Y. [Eep(^^')] [Y.e,{hy^)\e,{-hA{n + a)) 

^ 6=0 x=l y=0 

= P-^ + l E(^V^-0)^VP^p(-^^(^ + «)) (7.15) 

^ b=l P 

with C = 1 (resp. ( = i) ii p = l(mod 4) (resp. p = 3(mod 4)). 
Assume 

{n + a,p) = l. (7.16) 

For p = l(mod 4), 

p-i 

For p = 3 (mod 4), 

p-i 



1 U 

(7.15) = P - 2 - — E ~ ^^^^ + = P - 2 + a, |a| = 1. 



i h 

(7.15) =P - [-)ep{-bA{n + a)) ^ p + a' , \a'\ = 1. 

V ^ fe=l ^ 
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Hence, if (7.16), 

( p''(l--) if p = 1 (mod 4) 

(7.14) >\ ) 'I (7.17) 

1 pMl- M if p = 3 (mod 4). 

Consequently 

(''•12) ^^Yl l[(n+a,p) = l for 2<p<Qoo] H ("^ ~ ») 
J^' 3<p<Qoo 

and 

iy{n) > (logQoo)"^i?"^ l[(n+o,p) = l for 2<p<Qoo] (^-IS) 

where ^* refers to an additional congruence condition on n, f (mod 8). 

As an immediate consequence of (7.7) and (7.18), we obtain an alternative proof of 
the 'positive density' conjecture, first established in [B-Fl]. 

Theorem 1. Any integral Apollonian circle packing produces a set of curvatures of 
positive density in Z. 

Proof. 

With previous notation, let S C [1, P^R] n Z satisfy \S\ > (1 - t)P^R with r > a 
sufiiciently small constant. We show that 

J^a;(n)>0. (7.19) 

Since S is arbitrary, this will imply that [1,P^R] contains at least tP^R curvatures 
and hence Theorem 1. 



Let Qoo be a large constant and K — n2<p<Qoo^'" i'^-^^) 

1/2 

!y(n) - \S\''^[2^\u{n)-u{n)\^\ 

neS neS 



> (logQoo)-^ R~^J2Y1 M{n+a,K) = l] - P^Q 00 

T neS (7.20) 



(where (*) refers to a congruence condition (mod 8) on n + a). 
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Next 



neS n<P'^R 

>-\ n {^-ly^-rp'ii 

2<p<Qoo 

>{logQoo)-' P^R (7.21) 
for T = t(Qoo) small enough. Substituting in (7.20) gives 

J2 ^(n) > P^logQoo)-'' - Qoo^'°"'] > 

neS 

for QoQ large enough. 

Remark. It should be noted that the previous argument establishes a stronger state- 
ment in fact. It follows indeed that the set of curvatures in the ACP contains a coset 
of the integers, up to a 'small' exceptional set (where 'small' refers to small density). 
We do not attempt here to make quantitatively stronger statements, as better results 
will be obtained in the forthcoming paper [B-K2] . 

(8). Prime Curvatures 

Another application of our analysis is an analogue of Theorem 1 for primes. Thus 

Theorem 2. The set of curvatures produced by any primitive integral Apollonian 
circle packing contains a subset of the primes of positive density. 

Denote V the set of primes and assume 5 C [1, P^R] fl V satisfies 

X A(n) > (1 - T)p2i?. (8.1) 

n€S 

A straightforward adjustment of the proof of Theorem 1 shows that it will suffice to 
establish a bound 

I J2 A(n)[a;(n) - u{n)]\ < Q^^P^ (8.2) 
nes 

(c > some fixed constant, Qoo is a sufficiently large constant and r = t((5oo) > O). 
Defining 

T(^) = X A(n)e(n^) (8.3) 

neS 
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we have to bound 





< 1 \S^-v\.\T\ 


Jt 


Jt 



(8.4) 



Invoking an additional ingredient, we will rely on certain distributional properties of 
exponential sums of the form (8.3), assuming jS*! ~ P^R. 

As a consequence of the majorant property for the set of primes (see [B], [Gr]), 
there is the distributional inequality 

mes [e e T; \T{e)\ > SP'^R] <e 5"^"^(P^i?)"^ for S > (logP)"^ (8.5) 

(and hence for all 5 > 0). 

Hence, if C T such that P'^R\Q,\ > 1, we have 



/ \T{9)fd9<.{\n\P^RyP^R. 



(8.6) 



Take Qo = (logP)i°'° and let Mo{q,b),Moo{q,b) be defined by (6.4), (7.5). 
Decompose (8.4) as follows 



1(8.4)1 < / \S^-i>\ \T\ 

'T\ U Mo{q,b) 
q<Qo 



+ E 



/ 



Qoo<Q<Qo U Mo{q,b) 



Q dyadic 



q^Q 

{h,q)=l 



(8.7) 



(8.8) 



+ V / \S^-v\\T\ 

JMo{q,b) 

qSQoo 
{b,q)=l 

Recall that |i>| < l-S^^I by (7.4'). 

Since certainly ||T||2 < (log P)^/RP, we have 

(8.7) < (log P)-^VRP 



I 



r\ U Mo{q,b) 



(8.9) 



(6-7) 1 ^ P -iin-'^ 

< (log P) ^ Vp.P. ^ ^ '° < P2 (log 
y R 



(8.10) 
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We estimate (8.8). Fix Qoo < Q < Qo {Q dyadic) and define 



M'iq,b) 



e- 



< 



(8.11) 



Then 



\T\ < 



U Moiq,b) 
qr^Q 



j \S^\.\T\+ E / |5c.|.|T| = (8.12) + (8.13). 



U M'{q,h) 

By (6.7) and (8.6) 



q~Q Mo{q.h)\M'{q,b) 
{b,q) = l 



(8.12) < -^Q-^io-^Tl ^ ^,^^ J2 « -^Q-^''-'q'^pVr<Q-'s''-'p' (8.14) 

g~Q * 



and, using (6.5) and (8.6) 



dip 



q^Q 

{b,q) = l 



< 



s,Q2<2«<Q2 



d(fi 



q~Q |<p| — 2« 
(M)=i «^ 



P2 



and from (8.14), (8.15) 



(8.8) < Q 



00 



(8.15) 



(8.16) 



Finally, 
(8.9) < 



U U l\o-^\<§^. 

9<Q00 (b,q) = l 



= (8.17) + (8.18) 
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From (7.6), (8.6), 

(8.17) « il^^ - Ph Qoo {P^Rf^^ < ^ Qoo 
and (8.18) is bounded as (8.13) 



ilO" 



(8.18) < P^Qoo 



Hence 



From (8.10), (8.16), (8.21), 



(8.9) < Q 



GO 



) 3 



(8.4)<Q-5^^ P\ 



(8.19) 

(8.20) 
(8.21) 

(8.22) 



which is the desired inequahty (8.2). This proves Theorem 2. 
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